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ABSTRACT 

Let {v e (x, t )} e> o be a family of approximate solutions for the nonlinear scalar conservation 
law u t + f x (u) = 0 with Co-initial data. Assume that {u £ (a:, t )} are Lip + - stable in the sense 
that they satisfy Oleinik’s E-entropy condition. We prove that if these approximate solutions 
are Lip'- consist e n t , i.e., if ||u e (-,0) - u(-, 0)|j Ll y {x) + \\v £ t + fx(v e )\\ Lij> '^ t) = O(e), then they 
converge to the entropy solution and the convergence rate estimate, ||u £ (-,t) — w( - j<)ll£,«y(*) = 
0(e), holds. Consequently, sharp L p and pointwise error estimates are derived. 

We demonstrate these convergence rate results in the context of entropy satisfying finite- 
difference and Glimm’s schemes. 


1 This research was supported in part by NASA Contract No. NAS1-18605 while the authors were in 
residence at ICASE, NASA Langley Research Center, Hampton, VA 23665. Additional support was provided 
by ONR Contract No. N00014-91-J-1076. 




1. INTRODUCTION 


We are concerned here with the convergence rate of approximate solutions for the nonlin- 
ear scalar conservation law, u t + f x {u) = 0 with Co-initial data. In this context we first recall 
Strang’s theorem which shows that the classical Lax-Richtmyer linear convergence theory 
applies for such nonlinear problem, as long as the underlying solution is sufficiently smooth 
e.g., [Ri-Mo, §5]. Since the solutions of the nonlinear conservation law develop spontaneous 
shock-discontinuities at a finite time, Strang’s result does not apply beyond this critical 
time. Indeed, the Fourier method as well as other L 2 - conservative schemes provide simple 
counterexamples of a consistent approximations which fail to converge (to the discontinuous 
entropy solution), despite their linearized Instability, e.g., [Ta3],[Ta4]. 

In this paper we extend the linear convergence theory into the weak regime. The extension 
is based on the usual two ingredients of stability and consistency. On the one hand, the 
counterexamples mentioned above show that one must strengthen the linearized Instability 
requirement. We assume that the approximate solutions are Tip + -stable in the sense that 
they satisfy a one-sided Lipschitz condition, in agreement with Oleinik’s E-condition for 
the entropy solution. On the other hand, the lack of smoothness requires to weaken the 
consistency requirement, which is measured here in the Tip'- (semi) norm. In §2 we prove for 
I/zp + -stable approximate solutions, that their Lzp'-convergence rate to the entropy solution is 
of the same order as their Lip'-consistency. The Tzp'-convergence rate is then converted into 
stronger L p convergence rate estimates. In particular, we recover the usual L 1 -convergence 
rate of order |, and we obtain new sharp pointwise error estimates which depend on the 
local smoothness of the entropy solution. 

In §3 we implement these error estimates for finite-difference approximations, using a 
finite-element representation which is interesting for its own sake. In §4 we apply these error 
estimates for the Glimm scheme. Other applications of the current framework, to spectral 
viscosity approximations and various viscosity regularizations, can be found in- [Ta5],[Sc- 
Ta],[Ta6]. 


2. APPROXIMATE SOLUTIONS 

We study approximate solutions of the scalar, genuinely nonlinear conservation law 

r\ r\ 

(2.1) —u{x,t) + —f(u(x,t)) = 0, f">a> 0, 

with compactly supported initial conditions prescribed at / = 0, 

(2.2) u(x,t = 0) = u 0 (x). 
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Let {u c (a:, t)} £> o be a family of approximate solutions of the conservation law (2.1), (2.2) 
in the following sense. 

DEFINITION. 

A. We say that {u E (a:,t)} £ >o are conservative solutions if 

(2.3) f v e (x,t)dx = f Uo(x)dx, t > 0. 

J X J X 

B. We say that {i> e (a;,f)} £ >o are Lip' - consistent with the conservation law (2.1), (2.2) if 
the following estimates are fulfilled? : 

(i) consistency with the initial conditions (2.2), 

(2.4a) ||u e (a;, 0) - u 0 (a;)|| Ll y < I< 0 ■ e 

(ii) consistency with the conservation law (2.1), 

(2.46) |K(x, t) + fx(v e (x, t))|U,y(*,[o.T]) < K T • e. 

We are interested in the convergence rate of the approximate solutions, u £ (:r,f), as their 
small parameter e j 0. This requires an appropriate stability definition for such approximate 
solutions. Recall that the entropy solution of the nonlinear conservation law (2.1), (2.2) 
satisfies the a priori estimate [Br - Os], [Ta5] 

(2 -5) W-.()lb P+ < ||uo||z ^ +ct( , <>0. 

The case ||uo||i,ij>+ = oo is included in (2.5), and it corresponds to the exact ~ t~ l decay rate 
of an initial rarefaction. 

DEFINITION. We say that {v s (x, t)} £> o are Lip + - stable if there exists a constant ft > 0 
(independent oft and e) such that the following estimate, analogous to (2.5), is fulfilled: 

(2.6) IKO.OIIk,.* < | W .,0)|| 1 z j >t+/ K- ‘ 2 °- 


Remarks. 


2 We let ||^[|iip, ||^|U< P + and ||^|| L< y denote respectively, ess sup XJ t y 

and su p ^*AV, ’ where = /sup? * <?■ 


4>(x)-4>(y) 

x-y 




4>(x)-4>(y) 

x-y 


+ 
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(i) The case of an initial rarefaction subject to the quadratic flux f(u) = |u 2 demonstrates 
that the a priori decay estimate of the exact entropy solution in (2.5) is sharp. A 
comparison of (2.6) with (2.5) shows that a necessary condition for the convergence of 
{u £ } e>0 is 

(2.7) 0 < /? < a, 

for otherwise, the decay rate of {u £ (-,f)} (and hence of its e — > 0 limit) would be faster 
than that of the exact entropy solution. 

(ii) The case /? > 0 in (2.6) corresponds to a strict Tzp + -stability in the sense that 
||u £ (q f)||i t p+ decays in time, in agreement with the decay of rarefactions indicated 
in (2.5). 

(iii) In general, any a priori bound 

(2.8) ||v £ (-,t)||jr / ,p+ < Constr < oo, 0 < t < T, 

is a sufficient stability condition for the convergence results discussed below. In par- 
ticular, we allow for /? — 0 in (2.6), as long as the approximate initial conditions are 
Xzp + -bounded. We remark that the restriction of £ip + -bounded initial data is indeed 
necessary for convergence, in view of the counterexample of Roe’s scheme discussed in 
§3. Unless stated otherwise, we therefore restrict our attention to the class of Lip + - 
bounded (i.e., rarefaction-free) initial conditions, where 

(2.9) Lq = max(||uo||L ! >+,||u £ (-,0)|| iip +) < oo. 

Finally, we remark that in case of strict Tfp + -stability, i.e., in case (2.6) holds with 
/? > 0, then one can remove this restriction of -bounded initial data and our 
convergence results can be extended to include general L~ c -initial conditions. The 
discussion of this case will be dealt elsewhere. 

We begin with the following theorem which is at the heart of matter. 

THEOREM 2.1. 

A. Let {u £ (a;, t)} e >o be a family of conservative, Lip + -stable approximate solutions of the 
conservation law (2.1), (2. 2), subject to the Lip + -bounded initial conditions (2.9). Then the 
following error estimate holds 

(2.10a) |K(.,T) - u(-,T) \\ Lip , < C T [||u c (-,0) - u 0 (-)||l, p ' + IK + f x (v')\\ Lip , {Xt[0tT]) } , 
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where n 

C T ~(l+(3LtT)\ v = ^j-> 1. 

B. In particular, if the family {v c (x,t )} e > o is also Lip 1 -consistent of order 0(e), i.e., 
(2.4a), (2.4b) hold, then v e (x,t) converges to the entropy solution u(x,t ) and the following 
convergence rate estimate holds 

(2.106) T) — u{-,T)\\ Lip , < Mt • £, Mt = {Ko + Kt){ 1 + /3L%T) V . 

PROOF. We proceed along the lines of [Ta5]. The difference, e e (x,t) = v s {x,t) - u(x,t), 
satisfies the error equation 

(2.11) ^ e e {x,t ) + ^~[a s (x,t)e e (x,t) ] = F e (x,t), 

where a e (x,t ) stands for the mean-value 

a e (x,t)=f a[fv E (x,t) + (1 - Z)u(x,t)]df, a(-) = /'(•), 

and F s (x,t) is the truncation error, 

F e {x,t) = v E (x,t) + f x (v E {x,t)). 

Given an arbitrary <p(x)eWo'°°, we let {^ £ (a:,t)}o<<<T denote the solution of the backward 
transport equation 

(2.12a) <f> £ t(x,t) + a e (x,t)(t>l(x,t) = 0, t < T, 

corresponding to the endvalues, <f>(x), prescribed at t = T, 

(2.126) (f> e (x,T) = <f>(x). 

Here, the following a priori estimate holds [Ta5, Theorem 2.2] 

(2.13) Wi-itfhip < exp(^ ||a £ (-,T)|| L , p+ dr)||(^(x)|| iip , 0 < t < T. 

The Iip + -stability of the entropy solution (2.5) and its approximate solutions in (2.6), pro- 
vide us with the one-sided Lipschitz upper-bound required on the right-hand side of (2.13): 

(2.14) ||a,(-,r)||„,* < + ||u(-,r)|| K „ + ] < 
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Equipped with (2.13), (2.14) we conclude 


(2.15a) 


m-A iu < 

< CSr||^(*)||i,„ 0<t<T, C T = (1 + 8L+T)\ 

and employing (2.12a) we also have 

||^ £ (av)||L.'p[o,T] < |a|comaxo<t<r \\4> £ (-,t)\\ L ip(x) < 

(2.15 5) 

< \a\ooC T \\<l>(x)\\Lip, |a|oo = max |/'|. 

Of course, (2.12) is just the adjoint problem of the error equation (2.11) which gives us 


(2.16) (e'(; T), «•)) = (e*(-, 0), </>’(; 0)) + (F'(x, t), 4?(x, 

Conservation implies that eg = f e c (x,0 )dx = 0 and by (2.15a) we find 

|(e‘(.,0W(-,0))| < ||«*(-,0)|M|^(.,0)|| iip < 

< (l + /3i„ + r)"||e'(-,0)|| i ,y||^(x)|| tjp ; 


(2.17a) 


similarly, conservation implies that F§ = /i,[ 0 i t] F*(x, t)dxdt = 0 and by (2. 15a), (2. 15b) we 
find 

\(F e (x,t), $ ((E, 0)i 2 (®,[0,r])| ^ 11^ ( x i ^)IU«p'(a;,[ 0 ,T])|| < ^ ^)||Ltp(i:,[ 0 ,T]) < 

(2.175) 

< (1 + |a|oo)C , T||7 ?£ (a:,f)||L«p'(®,[o,r])!|^(^)||Lip- 

The error estimate (2.10a) follows from the last two estimates together with (2.16). □ 

The Ttp'-convergence rate estimate (2.10b) can be extended to more familiar Wf*- 
convergence rate estimates. The rest of this section is devoted to three Corollaries which 
summarize these extensions. 

We begin by noting that the conservation and Tip + -stability of v e (-,f) imply that v s (-,T) 
- and consequently that e e (-,T), have bounded variation, 

||e g (~,T)H BV < Const 

Using this, one can extend Theorem 2.1 into a general W~ 1,p error estimate (consult [Ta5, 
Theorem 5.1]); namely, there exists a constant (depending on Mr and || e c (- , T 1 ) ||,bv) such 
that 

p+1 

(2.185) ||i» e (-, T) - «(■, T)\\ w -i* < Const T • e , 1< p < oo. 
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The case p = 1 would correspond to the Lip'-e rror estimate (2.10). 

Theorem 2.1 also enables us to estimate the X p -convergence rate of the Zzp + -stable ap- 
proximate solutions { v £ (x , i)} s >o- To this end we use a Cg-unit mass modifier ((x) to denote 
<j>s = <f>* |C( j)- By Theorem 2.1 we have 

(2.19) |e'(.,r)»^(.)|<M T |Mz,-. 

This, together with the straightforward estimate (see e.g. [Ta5,§3]) 

(2.20) |e*(.,T) . W-) - *(-)]| < ||e*(.,T)|| B v«||^||z,», 

imply that for any compactly supported (freL 00 we have 

(2.21) \J^{x,T)4>{x)dx\< [Mr|-M|e £ (.,T)||B^]||^||Loo. 

Choosing the free parameter 6 ~ \fe, (2.21) with truncated <j> = [e £ (-,T)] p_1 yields 

(2.22) ||u £ (-, T) - u(-, T)H L p < Const T ■ e 5 ?, 1 < p < oo. 

Summarizing (2.18b) and (2.22) we state 

COROLLARY 2.2. Let {v e (x,t)} £> 0 be a family of conservative, Lip' - consistent andLip + - 
stable approximate solutions of the conservation law (2.1), (2.2), with Lip + -bounded initial 
conditions (2.9). Then the following convergence rate estimates hold 

(2.23) ||u £ (-,T) -u(-,T)\\ w -,,p < Const T ■ , 1 < p < oo, s = 0, 1. 

The error estimate (2.23) with ( s,p ) = (0,1) yields L 1 convergence rate of order 
which is familiar from the setup of monotone difference approximations [Ku], [Sa]. Of course, 
uniform convergence (which corresponds to (s,p) = (0, oo)) fails in this case, due to the 
possible presence of shock discontinuities in the entropy solution u(-,t). Instead, one seeks 
pointwise convergence away from the singular support of u(-,t). To this end, we employ a 
Cq(— 1, l)-unit mass modifier of the form (s(x) = jC(f )• The error estimate (2.10) asserts 
that 

(2.24) |(»*(-, T) *(,)(*) - (»(;T) * MWI < M T f 2 
Moreover, if £(x) is chosen so u.at 

(2.25a) Jx k ((x)dx = 0 for k = 1,2, . . . ,p — 1, 
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then a straightforward error estimate based on Taylor’s expansion yields 
(2.254) !(«(-, T) * &)(*) - u(x, T)| < ^||C|| t . • 

p\ 

where |u( p )| ;oc measures the degree of local smoothness of 

QP 

|w^|ioc = T)||£,<» c (^ +SU pp^). 

The last two inequalities imply 

COROLLARY 2.3. Let {v e (x,t )} e> o be a family of conservative, Lip' -consistent and Lip + - 
stable approximate solutions of the conservation law (2.1), (2.2), with Lip + -bounded initial 
conditions (2.9). Then, for any p-order mollifier (s( x ) = |C(f ) satisfying (2.25a), the 
following convergence rate estimate holds 

(2.26) |(u g (-,T) * &)(*) ~ u(x,T ) | < Const T (l + ^j --) . S ~ eife. 

p\ 

Corollary 2.3 shows that by post-processing the approximate solutions u £ (-, t ), we are able 
to recover the pointwise values of u(x , t) with an error as close to e as the local smoothness 
of u(-,t) permits. A similar treatment enables the recovery of the derivatives of u(x,t ) as 
well, consult [Ta5, §4]. 

The particular case p = 1 in (2.26), deserves special attention. In this case, post- 
processing of the approximate solution with arbitrary Cp-unit mass mollifier C(x), gives 
us 

(2.27) |(u £ (-, T) * Cfi)(«) - u i x , T) | < Const • (1 + !«*(-, T)\i oc ) • 6 ~ eK 

We claim that the pointwise convergence rate of order indicated in (2.27) holds 

even without post-processing of the approximate solution. Indeed, let us consider now the 
difference 

v s (x,T)-(v c (-,T)*( s )(x) = f[v'(x,T)-v s {x -y,T)]( s (y)dy = 

J y 

= [ f ”‘( x ’ — ~ — ~ - T ) | • -tC(| )<fy- 

Jy L -y \ o 6 

By choosing a positive Cj-unit mass mollifier (( x ) supported on (—1,0) then, thanks to the 
L 2 p + -stability condition (2.6), the integrand on the right does not exceed Const -6, and hence 

(2.28a) v*(x, T ) — (v e (-, T ) * G)( x ) 5; Const • 5 . 
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Similarly, a different choice of a positive Co-unit mass modifier £(x) supported on (0, 1) leads 
to 


(2.286) v £ (x, T ) - (u £ (-, T) * ( s )(x) > Const • 6. 

The last two inequalities (with 6 ~ £ 3 ) together with (2.27) show that the approximate 
solution itself converges with an ©(e^-rate, as asserted. 

We summarize what we have shown by stating the following. 


COROLLARY 2.4. Let {v'(x,t)} s> o be a family of conservative, Lip'- consistent andLip + - 
stable approximate solutions of the conservation law (2.1), (2.2), with Lip + -bounded initial 
conditions (2.9). Then the following convergence rate estimate holds 


(2.29) |v c (a;, T*) - u(x,T ) | < Const x , T • £», Const*, r ~ 1 + \ u x{-,T)\ Loa ^_ e ^ x+e $y 


Remark . The above derivation of pointwise error estimates applies in more general situations. 
Consider, for example, a family of approximate solutions, {v € (x, t)} c >o which satisfies a 
standard L l (rather than Lip') error estimate 

(2.30) l(v‘(;T) - u(;T), *(■))! < Constj- • 

Then our previous arguments show how to post-process v e (-,T) in order to recover the 
pointwise values of the entropy solution, u(x, T) with an error as close to \Je as the local 
smoothness of u(-,T) permits. In particular, using (2.30) with a positive Co-unit mass 
modifier, ^(x) = |C(f ) we obtain 

(2.31) |(»'(-, T) » c«)(») - M ;T) * 00MI < Const j. • ^HCIk- 
Using this together with 

|(n(-, T) * £$)(x) — u(x, T)| < (5||£||/,i • ||u j; ( > , ^Ollz-g^s+suppC)? 

we find 

(2.32) |(u(-,T) * £«)(x) - u(x,T)| < Constr(l -f |u x (-,T)|/ oc )£4 , 6~e*. 

If the approximate solutions {v £ (x, t)} E>0 are also Tip + -stable, then we may augment (2.32) 
with (2.28) to conclude the pointwise error estimate 

(2.33) |t> £ (x, T ) - u(x, T)| < Const*, T • e * , Const*, T ~ 1 + |«*(-, ^y 
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3. FINITE DIFFERENCE APPROXIMATIONS 


We want to solve the conservation law (2.1) - (2.2) by difference approximations. To this 
end we use a grid (x v = vA x,t n = nAt) with a fixed mesh-ratio A = ^ = Const. The 
approximate solution at these grid points, u" = w(xi,,f n ), is determined by a conservative 
difference approximation which takes the following viscosity form, e.g., [Tal] 


(3.i) »"+' = v ; - -[/(*,) - /(«;.,)] + W + , - K) - QUW - <-.)]. -so. 


and is subject to Xip + -bounded initial conditions, 


(3.2) V ° = ~KxL u °(Q d Z’ L o = IK||l,>+ < oo. 

Xp -i 

Let u*(x, t ) be the piecewise linear interpolant of our grid solution, u e (x„, t n ) = u", depending 
on the small discretization parameter e = Ax j 0. It is given by 


(3.3) = E^OM), A"(x,i) = A,(*)A m ('), 

j,m 

where A j(x) and A m (t) denote the usual ‘hat’ functions, 

A j(x) = ^min(x - Xj_ i,x j+1 - x)+, A m {t) = ^min(f - t m_ V m+1 _ f) + . 

To study the convergence rate of v Ax (x,t) as Ax J. 0, we first have to verify the conser- 
vation and the Xip'-consistency of the difference approximation. To this end we proceed as 
follows. 

We first note that v Ax (x,t) are clearly conservative, for by the choice of the initial con- 
ditions in (3.2), 

J v Ax (x, t)dx = ^ + K+i = + u °+i = / M x )dx. 

Moreover, these initial conditions are Lip 1 - consistent - in fact the following estimate which 
is left to the reader holds, (u A:c (x,0) — Uq(x) , (f>(x)) < Const • (Ax) 2 ||u 0 (^)||bv • ||<Kz)||l.> 
Finally, we turn to consider the Xip'-consistency with the conservation law (2.1). To this 
end we compare u A (x,t) with certain entropy conservative schemes constructed in [Ta2]. 

A straightforward computation (carried out in the Appendix) shows that there exists a 
bounded piecewise-constant function, D n (x) = 52j Dj + iXj+±( x )i (Xj+i-( x ) = characteristic 
function of (xj, Xj+i)), such that the difference approximation (3.1) recast into the equivalent 
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form 3 


d 


^-v Ax (x,t),K n u (x,t)\ + ^-f(v Ax (x,t)),A2(x,t) 


d 


,Ar/ 


dt 


(3.4) 


A x,t 


dx‘ 


x,A t 


Hence, for arbitrary <f>eC t q°, we may rewrite (3.4) as 


D(x),At 


(3.5) 




ETt- 

k = i 


The sum on the right-hand side of (3.5) represents the truncation error of the difference 
approximation (3.1), and according to (3.4), it consists of the following four contributions 

A 

(here, (f>(x,t) = J2t/, n ( f ) ( x i'iC)A^(x,t) denotes the piecewise-linear interpolant of cf>(x,t )): 

Tl = 

li = yluf'.AU 

T : , = (t>f',<S)„-(t-f*,#W, 

T< = A,. 


We want to show that the difference approximation (3.1) is consistent with the conservation 
law (2.1), in the sense that the Lip ’- size of its truncation error is of order 0( Ax). The 
required estimates in this direction are collected below. We begin with a straightforward 
estimate of the first term, 


(3.6a) 


\ T i\ ^ r |U 1 (|£>(*)I.A«) * ||<MU°°(*,At) 

< Ci ■ Ax||u Al (a:,t)|| L i( [0>r ] i sK( ; r)) • ||^(^,OIU*p(*,[o,T])- 


The difference approximation (3.1) enables us to upper bound time-differences in terms of 
spatial differences to yield the following upper-bound on the second term, 


(3.66) 


|r 2 1 < * II^IIl 00 ^) 

< C 2 • Aa:||u Ai: (a;,f)|| i i( [ 0 iT ] iB y( ;r )) • ||</>(M)ilL,>(^[o,T])- 


3 The Euclidean and weighted £ 2 -inner products are denoted by (p, V’la? = / p(x)rj;(x)dx and (p,iP)d(x) = 
J p(x)ip(x)D(x)dx. The corresponding discrete (?- inner product reads {p,ip) Ax = p{x v )4>{x v )Ax. 
Similar notations are used for (x, <)-functions, e.g., (p,iP)d^x), At = J2n f x P( x >t n )‘ l P( x >t n )D( x )dxAt, 
INIr, P( Ar,() = It \p(x„,t)\PAxdt, etc. 
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The third contribution to the truncation error we rewrite as 

T3 = [(n A:r , <i>) x ,t ~ (vf x , + ( u f <t> ~ 4 >)x,t = T31 + T 32 . 

We have (abbreviating )): 

v,j,n Z u,j,n Z 


Eh< +, -»r 1 ' 1 

i/,n " 


)g(« +I - 2« + «-.)Ax, 


and hence T 31 is upper bounded by 


I Tail < Const. Ax||u t Al || L i(A riA<) • ||<KM)||l«p(*,[o,t])- 
This together with the standard interpolation error estimate 

1^32 1 < |h A lLl(*,<) • M - j>\\L°°(x,t) < Const • AsHuf'lljr,!^) • ||^(x,t)||i,p( Xi [ 0 l T]), 
give us that the third term does not exceed 

|T 3 | < Const • Ax||v t Ax || L i( XtJ )||^(x,t)|| Lip ( Xi [ 0)?1 ) 

(3.6c) 


< C Z ■ Ax||u A *(x,f)|| L l ([ o iT ],BV(*)) ■ |MM)l|wp(*,[0,T])- 
A similar treatment of the fourth term implies 

(3.6d) \T 4 \ < C 4 • Ax||u Aa; (x,t)|| L i([o,T],BV(x)) • \\<f>(x,t)\\ L ip(x,[o,T])- 

Equipped with the last four estimates (3.6a) - (3.6d), we return to (3.5), obtaining 
(3.7) |(vf* + fx(v Ax ),(f>) x | < Const • Ax||u Aa; (x,t)|| L i( [ o t T],BV(x)) • \\<t>{x,t)\\ L ip(x,[o,T[y 


This shows that the Tip'-consistency estimate (2.4b) holds with e = Ax and 

K t ~ ||u A:r (x, <)||l 1 ([o,T],BV(*))* Thus, Corollaries 2.2 - 2.4 apply and their various error 

estimates are put together in the following. 


THEOREM 3.1. Assume that the difference approximation (3.1) -(3.2) is Lip + -stable in 
the sense that the following one-sided Lipschitz condition is fulfilled: 


(3.8) 


~2 < 


1 


Ax ~ [Lt]~ x + pt n ' 
Then the following error estimates hold: 


, 0 < t n < T, Au" + i = v v+x (t n ) - v v (t n ). 


(3.9a) 


u Ax (-, T) — u(-, T)||vy-«,p < Constj • (Ax) 2 p , 1 < p < 00 , s = 0, 1 


»p+i 
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v Ax (x, T ) — u(x, T)| < Const^r • (Ax) 3 , 


(3.96) 

where Const,. r ~ 1 + |«i(-, T)| , , 1 1 . 

X ' 1 x ' ’ yl Z,°°(a;-(Ai:)J ) 2r+(Aa;)S) 

EXAMPLES. The following first order accurate schemes (identified in an increasing order 
according to their numerical viscosity coefficient, Q„+i = Q" + i) are frequently referred to 
in the literature. 


(3.10a) Roe scheme : 


s\R \| | __ n f( v u- 

Qj/+§ — a u+\ — °W§ ~ A«” + j 


(3.106) Godunov scheme : 


Q° +l = Amax v 


/KUl)+/W)~ 2/(v) 

Av" , 

"+i 


(3.10c) Engquist - Osher scheme : Q &° , = ^ — f^ +1 \f(v)\dv 

2 

(3.10d) Lax — Friedrichs scheme : Q Lx f = 1. 

U '”2 

We comment briefly on the Lzp + -stability condition of these schemes. 

For the Roe (or Courant-Isaacson-Rees) scheme, Lip + -stability (3.8) with /? = 0 (no 
decay), was proved in [Br]. Note that the assumption of Lzp + -bounded initial conditions is 
essential for convergence to the entropy solution in this case, in view of the discrete steady- 
state solution, = sgn(^ + |), which shows that convergence of Roe scheme to the correct 
entropy rarefaction fails due to the Lzp + -unboundedness of the initial data. 

The Godunov and Lax-Friedrichs schemes can be viewed as cell averaging of the exact 
Riemann solver associated with (2.1), for which (2.5) holds: 


(3.11) 


L " +1 ~ [L+]- 1 + aAt' L *~ 


Au n i 

V+l 


Arguing along the lines of [Br - Os] we conclude from (3.11) that both Godunov and Lax- 
Friedrichs schemes satisfy the Lzp + -stability (and in fact the Lzp + -decay) (3.8) with /? = a 
in agreement with (2.8). One then recovers the convergence rate estimates quoted in the 
previous section, with error coefficients depending on Mj ~ (1 -f ccLq T) v , rj = . 

Finally, let us consider the Lzp + -stability of the Engquist-Osher scheme. It coincides with 
Godunov’s scheme except for sonic shock cells (where a(zv +1 ) < 0 < a(v v )), which leads to 


(3.12) 


A k+i| < Q»+ 1 < Qf+i < Q°+\ - Const • (Azv + i). 
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Hence the forward differences of E-0 scheme are upper-bounded by 

= (1 - + i(0«j - Aa„ +t )A< +} 




(3.13) 


+ - 0? + »)Av” +i + ;(«“* - CiKi 

-<5 

+ o(Q?-L + Aa„_i)Av"_i. 


1 

2 2 1 


< (1 - 0f + °|) A< +| + 1(0° } - Aa„ +| )A< +|+ 


1 

2 1 


We distinguish between two cases. If Au n i > 0, then the first term on the right of (3.13) 
does not exceed (1 — Q^ + i_)Au" +i , and hence the E-0 satisfies the one-sided Lipschitz 
condition in this case because Godunov’s scheme does. Otherwise, Au", i and therefore 
(1 — )Av* + i is negative, hence 

A y "+1 < i(0? + | - K + |)A< +} + 1 (<?° j 

and the Iip + -bound follows in view of (3.12) and the CFL condition A max |/'| < |. 

Using Theorem 3.1 we conclude 


COROLLARY 3.2. Consider the conservation law (2.1), (2.2) with Lip + -bounded initial 
data (2.9). Then the Roe, Godunov, Engquist-Osher, and Lax-Friedrichs difference approx- 
imations (3.1), (3.10) with discrete initial data (3.2) converge, and their piecewise-linear 
interpolants v Ax (x,t), satisfy the convergence rate estimates (3.9a), (3.9b). 


4. GLIMM SCHEME 

We recall the construction of Glimm approximate solution for the conservation law (2.1), 
see [Gl], [Sm]. We let v(x,t) be the entropy solution of (2.1) in the slabs t n <t < t n+1 ,n > 0. 
To proceed in time, the solution is extended with a jump discontinuity across the lines 
f n + x , n > 0, where v(x,f n+1 ) takes the piecewise constant values 

(4.1) v(x , r +1 ) = ^2 v(x v + 9 n Ax,t n+1 - 0)x I /(^)- 

V 

Notice that in each slab, v(x,t ) consists of successive noninteracting Riemann solutions 
provided the CFL condition, A • max |a(u)| < | is met. This defines the Glimm approximate 
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solution, v(x, t ) = i> £ (x, t), depending on the mesh parameters e = Ax = AA t, and the set of 
random variables {0 n }, uniformly distributed in [— |]. In the deterministic version of the 

Glimm scheme, Liu [Li] employs equidistributed rather than random sequence of numbers 
{0 n }. We note that in both cases, one makes use of exactly one random or equidistributed 
choice per time step (independently of the spatial cells), as was first advocated by Chorin 
[Cho]. This implies that both versions of Glimm scheme share the exact Lip + -decay of the 
entropy solution, for by (2.5) 

IM-> ( ” +, )IIkp* < - 0)[| tip+ < 

(4 ' 2) < 1 < L . 

+ “Af IM-,0)IU + + at"*' 

Namely, the Lzp + -stability (2.6) holds with f3 = a. 

Although Glimm approximate solutions are conservative “on the average,” they do not 
satisfy the conservation requirement (2.3). We therefore need to slightly modify our previous 
convergence arguments in this case. 


We first recall the truncation error estimate for the deterministic version of Glimm scheme 
[Ho - Sm, Theorem 3.2], 


(4.3) 


+ Mv Ax (x,t)), < f > ( x ^)) L2{Xt[0iT]) ^ 

< Const T [(Ax )2 [ log Ax| • ||</>||l~ + Ax • ||<^(a:,t)||i,ip( ;ri [o,r])] 


Let <f>(x,t ) = (f) Ax (x,t ) denote the solution of the adjoint error equation (2.12). Applying 
(4.3) instead of (2.17b) and arguing along the lines of Theorem 2.1, we conclude that Glimm 
scheme is Lip'-consistent (and hence has a Lip'-convergence rate) of order (Ax) 5 1 log Ax |, 


(4.4) |(e A;c (-, T), <?i(-))| < Const T [(Ax)* | log Ax| • |M|l~ + Ax • ||^(x)|| £ ,,-p] . 

To obtain an improved L p -convergence rate estimate we employ (4.4) with <f>s = (j) 5 

obtaining 


(4.5) 


l Ax ' 

|(e Ar (-,T),^)| < Const T (Ax) 5 1 log Ax| + — ||^(x)|| L oo 


Using this estimate (instead of (2.19)) together with (2.20) and choosing the free param- 
eter 8 ~ (Ax)*, we end up with an (almost) L p - convergence rate of order 0(Ax| log Axl) 5 ?. 
As noted before, the Lzp + -stability of Glimm’s approximate solutions enables us to convert 
the L 1 into pointwise convergence rate estimate. 
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We close this section by stating the following. 

THEOREM 4.1. Consider the conservation law (2.1), (2.2) with sufficiently small Lip + - 
bounded initial data (2.9). Then the (deterministic version of) Glimm approximate solution 
v Ax (x,t ) in (4.1) converges to the entropy solution u(x,t), and the following convergence rate 
estimates hold: 

(4.6) ||n A:r (-,T) — u(-,T)|| iP < Constr • (Ax|log Ai|)2p, 1 < p < oo, 

(4.7) |u Ax (x,T) — u(x,T) | < Const*,?’ • (Ax| log Ax|)*, 
where Const,, r ~ 1 + 

Remarks . 

1. A sharp Z^-error estimate of order 0(Ax)2 can be found in [Lu], improving the 
previous error estimates of [Ho-Sm]. 

2. Theorem 4.1 hinges on the truncation error estimate (4.3) which assumes initial data 
which sufficiently small variation [Ho - Sm]. Extensions to strong initial discontinuities for 
Glimm scheme and the front tracking method can be found in [Che, Theorems 4.6 and 5.2]. 
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APPENDIX 


We want to show that the piecewise-linear interpolant v Ax (x,t) in (3.3) serves as an 
approximate weak solution of the conservation law (2.1). 

Let 

v n ( x ) = ^u"Ai,(a;) and v„(t) = ]Pu£A n (t) 

V n 

denote the spatial and temporal interpolants of the discrete grid solution 
Straightforward integration by parts yields [Ta2] 

C U^ x {x,t)),K v (x)) x = |[/K + i(<)) -/K-! (<))]- 

(ia) 

— 2^i/+i ~ Q* v - 

where ( we abbreviate u„ + i(£,i) = \[v„{t) + u„ +1 (t)] + ^Av^t)) 

(“) = At,„ +I (t) • /“ d-n r(K +i ((,m. 

2 ^ 

In particular, for f(v) = v we have Q* = 0 and (la) yields 


(v£ x (x,t), A u (x)) x = J[u„ +1 (f) - w„_i(i)]. 


Exchanging the role of the x and t variables in the last equality we get 


( 2 ) 


t ), A ”(()), = i[«” +1 (x) - „-■(*)]. 


Moreover, with D(x) = T, 1/ D 1 /+ lx 1/+ i(x) we have 


( 3 ) 


Ax 


2 ( v x ( x >t)> (Ais(x))x)d(x) — -^[Dy+^Avy+iit) — D v _±Av v _i.(t)\ 


and by exchanging the role of the x and t variables in (3) we get 

« T {v ‘ x{x ' ‘)’( An (‘»‘>‘ = + »"-■(*)]. 

The equalities (1) - (4) imply 


(10 


= ^IfK+i) - /(Cl)] - ^[Q: + i(i n )A^ n + , - Q;_,(r)Ai;;u], 


(20 


Ax , 


(» r (*, 0 , a?(M))a*,< = =Rw n+1 - vr 1 ], 
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(3') 


Ax 


9 ( v x (h v ( x ,t)) x )D(x),&t — 9 [-^i/+i.Au", I D u _ lAu" l], 


At, 


(40 


At 


Ax, 


~(v?*(x,t), (A:( X ,t)) t ) Ax , t = -~K +i - 2v: - v:- 1 ]. 


The difference approximation (3.1) reads 

/\ / /\ n* 

(5) AxK +1 - <] = — j -[/«+ j ) - /(Cl)] + AC t ]. 

By (20 and (40, the left-hand side (LHS) of (5) equals 

LHS = ^K +1 -»r'] + ^[< +1 - 2< + c 1 ] = 

= (•’i 4 “’(i><),A;(x,())a I ,i- ^(t) t Al (x,(),(A;(x.i)) ( ) Ali< . 

Next, we set D n L = \Q n , ± — Q* , i{t n )] then by (N) and (30 the right-hand side (RHS) of 

l-'T 2 A 2 2 

(5) equals 

rhs = -y[/ft,-/K.,)] + j ( e^iA»; tr e;.iA'-; H ]+ 

K h a< + i - CjACi] = 


Ax 


= ~(/r(w (®,f)),A"(x,f)) 3;)At - — (u“(x,t),(A"(x,f)) x ) D(a;))At 


and (3.4) now follows. 
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